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ABSTRACT 


A continuum theory is presented for representing the thermoelastic behavior of 
composites that can be idealized as transversely isotropic. This theory is consistent with 
anisotropic viscoplastic theories being developed presently at NASA Lewis Research 
Center. A multiaxial statement of the theory is presented, as well as plane stress and plane 
strain reductions. Experimental determination of the required material parameters and 
their theoretical constraints are discussed. Simple homogeneously stressed elements are 
examined to illustrate the effect of fiber orientation on the resulting strain distribution. 
Finally, the multiaxial stress — strain relations are expressed in matrix form to simplify 
and accelerate implementation of the theory into structural analysis codes. 


INTRODUCTION 

Historically, the study of the mechanical behavior of fiber reinforced composite 
materials has been approached from two viewpoints, the microscopic and macroscopic. 

The microscopic view considers the constituents (fiber/matrix) separately, addressing in 
detail the interaction between individual fibers and surrounding matrix and the behavior at 
the fiber - matrix interface. Alternatively, the macroscopic approach considers the 
composite to be a material in its own right, (continuum), with its own experimentally 
measurable properties specified for the composite as a whole. This continuum approach 
provides a relatively efficient framework for the prediction of observed macroscopic 
deformation behavfor. Extensive research into both the micro and macroscopic viewpoints 
has been undertaken (e.g.[l,2]). 

Here, a continuum theory is presented for representing the thermoelastic behavior of 
composites that can be idealized as a homogeneous continua with locally definable 
directional characteristics. Although it is presumed here that a single preferential (fiber) 
direction is identifiable at each material point, thus admitting the idealization of local 
transverse isotropy, the theory is extendible to account for two (or more) identifiable 
preferential directions. 


As indicated in [3], homogenization of textured materials (e,g composites) and 
applicability of continuum mechanics depends relatively upon characteristic structural 
dimensions, the severity of gradients (i.e, stress, temperature, etc.), and the relative scale 
and periodicity of the internal structure of the material. Examination of these conditions 
reveals that for many anticipated aerospace applications of composites, the formulation of 
continuum based theories is justified. 

The objective of this study is the development of a transversely-isotropic 
thermoelastic theory which is consistent with present, continuum-based, anisotropic 
viscOplastic theories. Therefore the present work follows closely the earlier work of Spencer 
[2,4] and Robinson [3,5], and relies heavily upon invariant theory and the existence of a 
strain and complementary energy potential in the formulation of the constitutive 
equations. Major contributions of this study are 1) the recasting of earlier work in terms of 
physically meaningful total stress invariants, 2) discussion of a transversely isotropic 
multiaxial thermal strain tensor, 3) identification of a correspondence between plane stress 
and plane strain for a transversely isotropic material and 4) specification of an 
experimental program for the complete determination of the required elastic constants as 
well as a discussion of the theoretical restrictions on these parameters, 

This study begins with the multiaxial statement of the theory, followed by the plane 
stress and plane strain simplifications. Experimental determination of the required 
material parameters and their theoretical constraints are then discussed. Simple 
applications are then employed to illustrate the capability of the theory in representing 
transversely isotropic thermoelastic behavior. 


MULTIAXIAL STATEMENT OF THE THEORY 

Here constitutive equations are given for a linear hyperelastic solid reinforced by a 
single family of fibers, i.e. a locally transversely isotropic material. Given a hyperelastic 
material the stress and strain components are related through a normality structure 
utilizing either a strain energy or complementary energy function [2,6], i.e. 
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where 

W = Wtfy) (3) 

and 

0 « il(cr.j) (4; 

Here tr.j denotes the components of (Cauchy) stress, e- the components of infinitesimal 
mechanical strain. 

Transversely isotropic material symmetry is included in the potentials of equations 
(3) and (4) by introducing a directional tensor d-d., for example 4 . 

J 

W = W( £ij ,<i i d.) 
and 

n = tya.., d.d.) 

where now W and fi depend not only upon e-j and cr. . , respectively at a point but also 
upon the local fiber direction. The symmetric tensor d-dj is formed by a self product of the 
unit vector d- denoting the local fiber direction, as the sense of d- is immaterial. 

As W and fl each depend on two symmetric second order tensors, form invariance 
(objectivity) requires that they depend only on certain invariants and invariant products of 
their respective tensorial arguments (i.e., an integrity basis). The integrity basis, for a 
function comprised of two symmetric second order tensors, has ten invariants [7]. A subset 
consisting of four invariants (Ip^I^Ig) of the irreducible set of invariants is used, as dj is 
a unit vector and we desire W and U only to be quadratic functions (so the resulting stress 
strain relations are linear). Assuming a polynomial representation, the potential functions 
may be expressed as follows: 

W(e.j,d i d.) = A’PJ + B’P’ + C’P£ + D’P| + E’P£ (5) 
and 

n(^,d.dj) = A P x + B P 2 + C P 3 + D P 4 + E P 5 (6) 
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where 


p i =I i 


P 2 = r 2 “ r 5 + W l l -m! + «iU) 


P 3 = I 5- I 4 


( 7 ) 


p 4 = I 4 


and 


P 5 = V4 


I x = ,../3 


T.i ss 1 /2 o* • c- 
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L = d-d* cr- 
4 l J Ji 


( 8 ) 


h = d i d j ff jk ff ki 

An analogous set of invariants to those in equations (7) and ( 8 ) is chosen for W, i.e. 
PjJ,P£,...,P£ and by replacing tf- by in equations (7) and ( 8 ), 

The subset of total stress invariants employed in (7) are similar to those utilized by 
Robinson and Duffy [3] for the corresponding deviatoric invariants. These invariants 
correspond physically to ; P^ — the square of the mean (hydrostatic) stress state, P 2 - the 
square of the maximum transverse shear stress , Pg — the square of the maximum 
longitudinal shear stress , P^ — the square of the normal stress in the local fiber direction, 
and Pg — the product of the mean and normal stress in the local fiber direction. Figure 1 
defines the coordinate system and P 2 ,Pg <md P 4 schematically. 



Substituting equations (5) - (8) into equations (3) and (4) result in the following 
linear elastic stress strain relations', 

ffj, = l/6{ (4A’~ 9B’) I{ + (3B’ + 2E’)I| } f. . n 

+ 1/2{(3B’ + 2E’)I j + [B’+4(D’~C*)]*4 } djdj 

+ B’Cjj + (C’ — B’)aJ, (9) 

and 

e-j = l/6{ (4A- 9B) I, + (3B + 2E)I 4 } . 

+ 1/2{(3B + 2E)Ij + [B + 4(D -C))I 4 } djdj 

+ Bo ij + (C-B)a ij (10) 


where 

6 - • - denotes the Kronecker delta function 

a ij “ ^i^jk^k + ^k^ki^j 
a ij = d i € jk d k fd ^ki d j 

Clearly equations (9) and (10) posses identical form; the coefficients and invariants, 
he wever are distinctly different. The unprimed coefficients A,B,C,D and E and their 
primed counterparts may be associated with more physical parameters (e.g. Young’s 
modulus and Poisson ratio) by conducting various thought experiments, For example 
consider the four stress states depicted in figure 2, given a preferred direction dj along the 2 
axis, i.e. d=(0,l,0). 


iThe appendix contains identical expressions to those of (9) and (10), yet written in matrix 
notation for easy implementation into structural analysis codes. 



From these states of stress (tests), it is easily shown that 


test a) 




B - *13/013 - 1 / 2 G t 

(11) 

test b) 


0 = h2 /a 12 = 1/2G l 

( 12 ) 

test c) 


2/9(A + 9D + 3E) = e 22 /a 22 « l/E L 

(13) 


- 1/9(2 A + 3E)E l » - e, ,/e 22 = v L 

(14) 


and 

test d) 


2(A/9 + B/4) = e, j/nr, , = 1/E T 

(15) 

-(4A - 9B) / f|| = " T 

(16) 

4 A + 9B 


-2(A/9 + E/ 6 )E t = - 633 / 6 ,! = ul 

(17) 


Where 

Grj, — Shear modulus (transverse) for the plane of isotropy 

Gj - Shear modulus (longitudinal) for a plane normal to the plane of isotropy, 

E«j» — Young’s modulus (transverse) in the plane of isotropy. 

E^~ Young’s modulus (longitudinal) normal to the plane of isotropy. 

Vy— Poisson’s ratio (transverse) that characterizes the transverse strain reduction 
in the plane of isotropy due to a tensile stress in the same plane. 

Ujj- Poisson’s ratio (longitudinal) that characterizes the transverse strain 

reduction in the plane of isotropy due to a tensile stress in a direction normal 
to it. 
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vfc- Poisson’s ratio that characterizes the strain reduction normal to the plane of 
isotropy due to a tensile stress in the plane of isotropy, 


Solving equations (13) thru (17) we find 



B = l/2Grp 
C = 1/2 G l 

D = (1+2*/ l )/2B Ij + (1-^ t )/4E t 


(18) 


E •= ~3{ (1— i/ t )/2E t + } 

and ; v ;| .. . 

2G«p ~ Eiji/(14'^<ji) 

v i = i, l e t/ e l 


Note that the five independent constants are chosen to be and G^. A 

process can be repeated to obtain the primed coefficients in terms of these five 
independent physical parameters, they are: 

A’ = 

B’ = 2G t 


U,E u /[ 


E L (i- l .r)-2E T ,J 



C- = 2G l (19) 

D' = E l {1 + E T (l-2i, L )V(2(E L (l-^)-2E T ^)]}/2 
E’ = -3 {E t E l (1-2i/ l )/[E l (1- 1 / t )- 2E tI /£]}/2 
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Therefore, the stress-strain relations of (9) and (10) may be expressed either in terms of 
the nonphysical coefficients A ( B,C,D and E, or alternatively in terms of the five 
independent physical parameters discussed above. 

Agreement between the present work and earlier work by Spencer [4] is evident, 
provided appropriate substitutions are made for the coefficients in equation (9) prior to 
comparison of the strain energy formulations. Additionally, a direct comparison with 
Lekhnitskii’s [8] work can be made when the physical parameters expressed in (18) are 
inserted into equation (10) and the preferred direction is assumed to be d— (0,0,1). 
Further, with appropriate substitution of equations (18) and (19) into (9) and (10), and 
considering conditions of isotropy, i.e. 


E = E L = E T ; v- 


T , 

d.d.~ &./3 
l j V 


T 


and 


( 20 ) 


the classical isotropic linear elastic stress — strain relations are obtained. 

The transversely isotropic stress-strain relations given above may be extended to 
include the effect of temperature by applying the traditional assumption of the additive 
nature of mechanical and thermal strain, i.e. 


where 





e- - total strain 
c.j - mechanical strain 


( 21 ) 


and the multiaxial thermal strain is assumed to have the following form 


<‘j = K“L-“T) d i d j + VtI AT (22) 

where a^~ characterizes the thermal expansion normal to the plane of isotropy 

Oj,— characterizes the thermal expansion in the plane of isotropy 

AT as T 2 - T x ; TpTg are two distinct temperatures 
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A similar expression to that of (22) was derived by Chadwick and Seet [9] from the 
Helmhotz free energy for a heat conducting elastic transversely isotropic material, There it 
was also shown that such a form satisfies thermodynamic restrictions, 

Again, imposing the conditions of isotropy (a=o^~a,p) it is evident that (22) 
reduces to 


ef.^&.aAT (23) 

the classical isotropic thermal strain tensor. Equation (23) is restricted to an isotropic 
continuum, as it implies that a change in temperature causes only a change in volume. This 
is not necessarily the case in (22), in that for an off axis fiber orientation a change in shape 
(shear) will occur with a change in temperature. An example of this is provided in a later 
section, 


TWO DIMENSIONAL SIMPLIFICATION 

Here, simplification of the multiaxial linear thermoelastic stress - strain relations of 
(21), in which equations (10) and (22) are incorporated, is made to that of plane stress and 
plane strain. For convenience equations (18), representing the relationship between 
physical and nonphysical parameters are substituted into (10). Plane stress, in the plane of 
isotropy (see figure 3), requires that 


T i = Vj 


0 


while plane strain requires 


3 i = e ij n j = 0 


where T- and 3. represent the stress and strain traction, respectively. If the local fiber 
direction (d.) and nj (indicating the directionality of the plane stress or strain assumption) 
are defined as, 


dj =* cos 0 

d 0 = sin 0 

4=o 


n, = 0 

n 2 “ 0 
n 3“ 1 
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the following expressions relating the nonzero stress and strain components are obtained. 


Plane Stress; 

e U ss<r u/ E T“' ,/ t <r 22 / E L + 

[l/2G Ii ~l/E ( y~^/E^)(2£r 11 co8 2 tl + 2<7 12 cos0sin0) + 

[l/E^+(l+ 2 ^)/E^"*l/Gj ( ](<r^co 8 2 0 + 

o^sin 2 0-f 2<r i2 cos &in 0)cos 2 0 

+ ( or»j. + (« L ~<^)cos 2 0|AT 


( 24 ) 


e 22 “ 


ff 22/ E T - Vll l\ + 

[l/2Gjv“l/E,j,~i/L/E L )(2£r 22 8in 2 0 4- 2<r 12 cos0sin0) + 
[i/E^ 4 , (l+ 2 ^)/Ej j “«l/G^](cr^^cO 8 2 0 4 
<r 22 sin 2 0 4 2o'j 2 co8&in0)sin 2 0 
4 [ a T 4 (a L -a T ) 8 in 2 0 JAT 


(25) 


©12 “ ^12 /^L "** IVSGl-I/Et-^/ElK cr n + O‘ 22 )cos0sm0 + 

[l/Erp+(l-f- 2 ^)/E^— l/Gj j ](cr^cos 2 0 4 
<r 22 sin 2 0+2 cr 12 cos 0sin0)cos 0sin 0 
4 [(a^--arji)co808in0jAT (26) 

e 33 “ ” v T/ E T^ll +<r 22) + 

[v t /E t ~u l IE l ]( <7 ^ ^cos 2 0-f ^ 22 sin 2 0+ 2 o- 12 cos ^in 0) 

4 a T AT (27) 
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Plane Strain; 


e ll = ( l ~ u i) a vJ E T ~ 22/% + 

[l/2Gr (l+t/mJ/Ei* 1(2^ . cos 2 0 4* 2(r. 9 cosflsin0) 

+ [(l“< / rji 9 )/B(p4 [1+2^(1+ , ]/^““V®£ j ]( <7 'ii cos2 ^ 

4* <r 22 8in 2 0+2<r 12 co80sin0)p,os 2 0 

4* {(1+^)^ 4- o&jjcos 2 0} AT (28) 

e 22 88 (l~ ll/ ^) <r 22/®T *” “** 

[l/2G^l“^)/B^-i/ L (l4*^ T )/E L ](2^ 22 sin 2 0 4* 2<r 12 cos0sin0) 

4- [(!“X / »ji 2 )/Eiji4'[14'2^(14'* / 'rji)“"^j ) * / j i , ]/Ej i p ”l/GjJ(<r^co8 2 0 
4* <r 22 sin 2 04*2<r 12 cos0sin0)8in 2 0 

4- {(l4-^ T )a T 4* [a I< ~(l4*^ T “'^£)0( T ]sin 2 d}AT (29) 

e 12 88 W 20 ! + 

[1/2Gl^1^-J)/E»ji*~i/jv(14-^t)/®lI( 0-^4* ^22) C08 ^* n ^ 
[(1 “^)/E ( j4’[14’2^(14'^(ji) ~ ^^£]/EjQ--i/Gj^3(c r 22 C082 ^ 

<r 22 sin 2 04*2 cr 1 2 cos 0sin0)cos O&inO 

+ { [ a^H! 1 4- ) OjJ eos flsin 0} A T (30) 

*33 = ^ T (^ u +<r 22 ) - [v T ~ vfj 

( <7j, jCOS 2 04- tr^sm 2 <?4-2o- 12 cos 0sin0) 

- a r E T AT (31) 

Clearly in order to convert an isothermal plane stress problem into an isothermal 
plane strain problem, provided nj is in the plane of isotropy (n.d. = 0), the following 
substitutions should be made: 



E x » E t / (t~v%) 


( 32 ) 


v jj 88 

where as before 



Employing the assumptions of isotropy, i.e, E^E^-E^ and the above reduce to 

their isotropic counterparts, as they must. Considering the nonisothermal case, the lack of 
a unique correspondence in the transverse coefficient of thermal expansion (o^) is noted 
when comparing the plane stress and plane strain expressions. This lack of correspondence, 
arises due to the difference in strain reduction, i.e. poisson ratio, for in plane and out of 
plane loadings. Clearly however, under isotropic conditions a functional correspondence for 
the thermal expansion coefficient does exist, i.e. a(l+i'). 

Alternatively, the assumptions of plane stress or strain may be imposed out of the 
plane of isotropy (n-d.^ 0). Here, in contrast to the previous assumptions taken in the 
plane of isotropy, the conversion from a plane stress to a plane strain problem is now a 
function of fiber orientation, or angle 0, Consider for example the case when n = (1,0,0) 
and 0-0° (i.e. the 1 axis is normal to the plane of isotropy) the following substitutions are 
required, 


Erp = E t /(1-^) 

"t = (i/ t + Vl>/( 1 -v£) 


(33) 


while if 0=90° (i.e. the 1 axis is in the plane of isotropy) identical substitutions as those in 
(32) are determined, as expected. For brevity, the angle dependent plane stress and plane 
strain expressions are omitted here. 
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EXPERIMENTAL DETERMINATION OF MATERIAL PARAMETERS 

Two types of specimens are presumed to be available. These consist of thin walled 
composite tubes that are longitudinally reinforced (having a single fiber direction oriented 
axially) and those that are helically reinforced, with a fiber orientation of 45<> (see figure 4). 
Each specimen may be loaded either in tension and/or torsion. 

The motivation for utilizing thin walled tubular specimens is two fold. The first is 
the fact that a thin walled tube is an ideal specimen for the development of constitutive 
relationships as it provides a nearly homogeneous region of stress and strain, and is 
statically determinate. Secondly, identical specimens may be employed, as well, to 
characterize inelastic material parameters [3,5] and combined tension/torsion experiments 
can be used as verification tests to assess the correctness of the multiaxial theory [ 10 ], both 
elastically and inelastically. 

First consider the longitudinal reinforced specimen, i.e. 0- 0°, subjected to an axial 
tensile load (o^). Clearly, then from equations (24) and (25) we obtain 


E L ~ *11^11 
and 

\ = ~ c 22^ c ll 


(34) 

(35) 


where is the applied axial stress and e jp «22 are measure d axial and circumferential 
strains respectively. Also subjecting the longitudinal specimen to a pure torsional load 
results in the quantification of the longitudinal shear modulus, as evident from equation 
( 26 ), 


G r. = »Wt 


12 ' '12 


(36) 


where 7 ^ is the measured engineering shear strain ( 2 e 12 ) and cr 12 th e applied torsional 
stress. 

Both the longitudinal and transverse thermal expansion coefficients also may be 
obtained by imposing a uniform temperature excursion, from some reference temperature, 
on the longitudinally reinforced specimen. Assuming no additional loads are applied to the 
specimen, it is easily seen from equations (24) and (25) that 

«!_ = €*, ,/AT 
and 

a T ” c 22/^ 
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where AT is known and e j j and are measured quantities, ; 

Now consider the helically reinforced specimen (Fig, 4), with a fiber orientation of 
450, subjected to a state of pure torsion. Equations (24) thru (27) then simplify to, 

“ e 22 ” - 1/2E^,)(7 1 2 (27) 


f 12 = (l/2Erj, (l+2*/k)/2EjJ<?^2 (28) 

and 

e 33 ~ (Mj/Et "** / L/®L) <r i2 ( 3 ®) 

from which the transverse modulus and poisson ratio of the composite material can be 

deduced. For example, from equations (37) and (39) we obtain 

E t » E L /(H2e u E L /cr 12 )] (40) 

and 

U T 88 ^ 3^12 + p L$l) ( 41 ) 


while equation (38) allows verification of either E^,^, or E^. Note that equation (41) 
assumes that the radial (diametral) strain (e^) is a measurable quantity. This 
measurement, however, is at best difficult. Clearly, if this measurement proves to elusive, 
an alternative test will need to be conducted to determine either the transverse shear 
modulus (Grji) or the transverse poisson ratio ( 14 -p). As an example, the double shear test 
may be employed to obtain a measure of the transverse shear modulus [11] . 

Although, the author prefers the use of tubular specimens, the high cost of 
fabrication and required test equipment may be prohibitive. Thus alternative test 
procedures, utilizing plate specimens may also be employed to obtain the required 
independent material parameters. The four material parameters E^,E'^, and can be 
obtained by performing tensile tests on plate specimens with orientations of 0 and 90°, see 
figure 5. Difficulties arise however, when attempting to find G^, as in~plane shear test 
methods must be employed; such as rail shear, symmetric rail shear, + 45° tension coupon, 
off axis tensile coupon, ARCAN and lOSIPESCU. The symmetric rail shear and + 45° 
tension coupon appear to be the preferred methods [12]. 

Finally, it is important to realize that although the above tests are viewed as 
characterization tests for the present theory, they also may be viewed as verification tests 
for the various micromechanics theories. 
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RESTRICTIONS ON THE MATERIAL PARAMETERS 

Theoretical constraints, which insure the admissibility of the elastic constants 
incorporated in the above expressions, 

may be derived by requiring that the complementary energy (strain energy) potential be 
positive definite; thereby guaranteeing convexity and uniqueness of solution. These 
inequalities are 

e l >0 

g l >0 

G t > 0 (42) 

V e t > v i 

1_I/ T > 2i '£ E t/ E L 

Under isotropic conditions (i.e. E=Ejy==E,p G=Gj^=G,j, and v=Vj—Vrj) these inequalities 
simplify to the well known restrictions imposed on elastic isotropic materials [6]. 

An additional restriction is imposed on the longitudinal shear modulus which 
insures that the global longitudinal (transverse) stiffness modulus be monotonically 
decreasing (increasing) for all fiber orientations between zero and ninety degrees. The 
necessary inequality is 


G l < E L /[2(l+v L )j 


(43) 


and may be deduced from either of the following conditions: 


dE n /d0=4e 12 <r n /e 2 n <0 

(44) 

dE 22 /d0 = "~ 4e l2 ^22^22 > 0 

(45) 


where ^n^22 are S^bal longitudinal and transverse stiffness moduli respectively and 

e ll ,€ 22 an< * e 12 are given e ^ uations (24) thru (26). 

Figures 6 and 7 illustrate the resulting variation in stiffness (for sample material 1 
and 2, respectively) along the 1 axis with respect to fiber orientation when the inequality 
of equation (43) is violated and satisfied, respectively. Table 1 provides the two complete 
sets of example material properties. Clearly, when G^ violates (43) a non realistic 
variation in stiffness, with respect to fiber orientation, results (see figure 5), The increase 


15 



and subsequent decrease in stiffness is attributed to the change in sign of the shear strain, 
as illustrated in figure 8, and implied in (44). In figure 8, the corresponding normal and 
shear strain variations verses fiber orientation are displayed, utilizing the material 
properties of sample one. Alternatively, figure 9 shows the normal and shear strain 
variations associated with the material parameters of sample two, i.e, figure 7. Also the 
importance of the sign, as well as magnitude, of the shear strain, i.e, c 12 , is evident from 
the expressions (44) and (45). 

The inequalities in equations (42) and (43) not only constrain the values of the 
elastic constants in the mathematical formulation to agree with certain basic physical 
principles, but also provide a means to determine whether experimentally obtained data is 
physically consistent with the mathematical model. If the measured material properties 
satisfy the constraints, one can proceed with confidence. Otherwise one might doubt the 
experimental techniques employed; or alternatively, the validity of the assumed 
idealization of transversely isotropic behavior for this particular material. 

SIMPLE ILLUSTRATIVE STATES OF STRESS 

Simple homogeneously stressed elements are employed in this section to illustrate 
the effect of fiber orientation on the resulting strain distributions. All results are 
associated with the material properties of sample 2 in Table 1. Four mechanical stress 
States and one thermal load are addressed. These are a uniaxial state, a pure shear state, 
an equal and opposite biaxial state, a two to one biaxial state, and a change in temperature 
with no applied stress. The resulting variations in strain with fiber orientation are shown 
in figures 9 thru 13, respectively. 

Examination of these figures reveals that the strain response , as one might expect, 
is highly dependent upon fiber orientation. Further, qualitative agreement with physical 
reasoning and experimentally observed response is also noted, thus providing additional 
confidence in the correctness of the present thermoelastic stress-strain relations. For 
example, consider the uniaxial stress state of figure 9. There the strain component is 
seen to be extensional (positive) for all fiber orientations, while e 22 is initially contractive 
(negative, e.g. at 0=0°) yet reaches a maximum extensional magnitude at 0=45°. This 
change in sign for the e 20 strain component would be completely unanticipated if one 
applied intuition developed from isotropic materials. Additionally no contraction or 
extension along the 2 axis would be observed when the fibers are orientated at either 0= 
17.05° or 72.95°, see figure 9. Such points of interest may be examined experimentally to 
verify the theory. 


Examining figures 10 and 11 (the pure shear and equal and opposite biaxial stress 
state, respectively) one must realize that although the strain responses appear to be 
distinct, they are in fact equivalent; as expected since the applied stress states are 
equivalent, provided the proper transformation is applied, For example with a rotation of 
45°, the biaxial stress state of figure 11, becomes one of pure shear, thus suggesting that the 
strain response be equivalent to that of figure 10 at 0 = 0 °, Indeed, taking into account a 
45® transformation, this can be easily shown to be the case. 

Consider the two to one biaxial stress state given in figure 12. This stress state is 
similar to that which arises in a thin walled tube with closed ends under internal pressure. 
As expected, the hoop strain (e 22 ) is significantly larger than the axial strain (c-q) at 0=O<> 
(a longitudinally reinforced tube). Increasing the fiber orientation toward 90° (i.e. 
circumferentially reinforced) decreases the hoop strain while increasing the axial. 
Interestingly, beyond 55° (a material dependent angle measurement) the magnitude of the 
axial strain surpasses that of the hoop strain, where at 90° is approximately twice that 
of e 2 2 > even though cr 22 = 2 cr^. One obvious application of the present theory is the 
optimization of fiber orientation, such that specific design criteria are satisfied. Adoption 
of known optimization principles to the present theory is straight forward: due to the 
inclusion of directional dependence (material symmetry) in the potential function. 

The effect of fiber orientation on the variation in thermal strains is illustrated in 
figure 13. Unlike the isotropic case, here a change in temperature may produce (for off 
axis fiber orientations) significant shear strains in addition to volumetric strains, The 
maximum shear strain, under a purely thermal load, will always occur at an angle of 45°, 
independent of material constants. 
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SUMMARY 

A continuum theory has been presented that represents the linear thermoelastic 
behavior of materials which can be idealized as transversely isotropic, Specifically this 
^theory may be applied to composite materials with a single preferred direction; wherein the 
composite is considered a material in its own right, with its own experimentally measurable 
properties specified for the composite as a whole. 

Some of the main points of the present study are: 


1) Tpe present theory is shown to be in agreement with earlier work done by 
Spencer, Lekhnitskii, and Chadwick and Seet. Also, assuming isotropic 
conditions, the present theory is shown to simplify to the classical 
thermoelastic stress-strain relations. 

2) Correspondence between isothermal plane stress and plane strain problems 
was established, while for the assumed multiaxial form of the thermal strain 
no correspondence is found for the nonisothermal case. 

3) Specification of an experimental procedure for the complete determination of 
the required material parameters is presented. Thin walled tubes, 
longitudinally and helically reinforced, loaded in tension and torsion are the 
primary specimens. 

4) Theoretical constraints on the values of the elastic constants, which insure 
their admissibility, are discussed. 

5} Simple illustrative states of stress indicate that the resulting strain response, 
in the presence of material anisotropy, may be highly nonintuitive. 
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APPENDIX - MATRIX REPRESENTATION OF MULTI AXIAL FORMULATION 

The multiaxial stress-strain relations given in index notation in equations (9) and 
(10) are expanded and expressed here in matrix notation. This is done in order to simplify 
and accelerate implementation of the present theory into structural analysis codes. 
Equations (9) and (10) now become 

M - [C] W (a) 

and 

W = [DJ M (b) 

respectively. Where the column matrices {a} and {e} contain the six independent tensorial 
components of o - and e... That is 

•j 

T 

M = l <T H> cf 22 ,cr 33 ,<T 12 ,Gr 23 ,<r 31^ 
and 

W “ t € ll’ e 22» e 33* € 12» € 23> e 3l] 

Note the inclusion of thermal strains, in the above equations, is accomplished by utilizing 
equation (22), i.e. substituting 


into equations (a) and (b). 

The independent components of the symmetric stiffness matrix, [C], are now 
expressed in terms of the nonphysical parameters as follows: 


* 
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First Row: 

c n = 2 A’/9 + B’/2 + [2C’-B’+2E’/3)dj + [B>/2 + 2(D'-C')]dJ 
C 12 = 2A’/9 - B'/2 + [B’/2 +E-/3](d|+dJ) + [B'/2 + 2(D’-C')]d jdj 
C 13 = 2A>/9 - B’/2 + [B'/2 +E>/3](dJ+dJ) + [B’/2 + 2(D’-C’))d}dJ 
C u = [C--B'/2+E'/3]d 1 d 2 + (B’/2 + 2(D'-C')]djd 2 
C ls = [B’/2+E73]d 2 d 3 + [B>/2 + 2(D’-C’))d}i,d 3 
C 16 = [C’-B’/2+E'/3]d 1 d 3 + [B’/2 + 2(D>-C’)]d|d 3 
Second Row: 

C 22 = 2A’/9 + B'/2 + [2C’-B’+2E’/3]dj + [B’/2 + 2(D’~C>)]d| 

C 23 = 2A'/9 - B’/2 + [B’/2 +E’/3](dJ+dJ) + [B’/2 + 2(D’-C>)Idgd§ 

C 24 = [0’-B>/2+E’/3]d 1 d 2 + [B>/2 + 2(D’-C’)]dJd 1 
C 2g = [C’-B'/2+E’/3]d 2 d 3 + [B’/2 + 2(D>-C>)]djd 3 
C 2g = [B'/2+E'/3)d 1 d 3 + [B>/2 + 2(D'-C’)]dJdgd 1 
Third Row: 

C 33 = 2A’/9 + B’/2 + [2C’-B’+2E’/3]d5 + [B'/2 + 2(D>-C’)]dJ 
C 34 = P 72+E73]d 1 d 2 + [B72 + 2(D’-C>)3d§d 2 d l 
Cg 5 = [C’-B72+E73]d 2 dg + [B72 + 2(D»-C*)]dJd 2 
Cg 6 = [B72+E73]d 1 dg + [B72 + 2(D’-C»)]d^d 1 
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Fourth Row: 


c 44 = 

B- + [O’— B’)(dj + dj) 

+ [B’/2 + 2(D’-C’)]djdj 

C 45 ~ 

[C’-B’Jdjdj + [B’/2 + 

2(D’-C’)]djdjdj 

C 46“ 

[C’-B’Jdjdj + [B’/2 + 

2(D’-C’)]d{d 2 d 3 

Fifth Row: 



C 55 = 

B’ + [C’-B’](dj + dj) 

+ [B>/2 + 2(D’-C’))djdjj 

C 56 = 

[C’-B’]d 1 d 2 + [B’/2 + 

2(D’-C’)]djd 2 dJ 

Sixth Row: 



C 66 = 

B’ + [O’— B’](dJ + dj) 

+ [B’/2 + 2(D’-C’)]d}d§ 


The components of the Compliance matrix, [D], are identical to those above, provided the 
unprimed coefficients A,B,C,D and E are substituted for their primed counterparts. For 
convenience the definitions of these nonphysical parameters, in terms of the more physical 
parameters, are repeated (c.f. equations (18) and (19)). 


The unprimed coefficients are 
A = 9/4[(l~^)/E^,] 

B = 1/2G t 

C = 1/2G l 

D = (1+2i/j^)/2Ej j -f (1— z/(j')/4E/j, 
E = — 3{ (1-j/ t )/2E t + } 

where 

2Gt = E/ji/(l*t"t / iji) 

v i = “tV®! 
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// , ■ 

and the primed are 

A> = 9/4 jB T E t /[E L (l-v T ) - 2E T t/£)| 

B' = 2G t 

C' = 2G l (d) ' 

D- = E l {1 + E t (1-2 1 / l )V|2(E l (1-./ t )- 2E t ^)]}/2 
E> = -3 {E t E I/ (1-2 1 / l )/[E l (1- 1 , t )- 2E t ^]}/2 

Table 1 

i n k s i * 1 0 -6 

material Ej^ Gj^ Mm ot^ Qbi 

sample X 28000 13000 18000 0.25 0.35 15.0 3.C 

sample 2 50000 13000 15000 0.25 0.35 15.0 3.0 





FIGURE 1, - DEFINITION OF COORDINATE SYSTEM AND 
PHYSICAL INVARIANTS KITH RESPECT TO THE PRE- 
FERRED FIBER DIRECTION d, 






FIGURE % - FUNDAMENTAL STRESS 
STATES, WITH RESPECT TO A GIVEN 
FIBER DIRECTION, REQUIRED TO 
DETERMINE MATERIAL PARAMETERS. 


r LONGITUDINAL 



FIGURE 3. - DEFINITION OF COORDINATE SYSTEM AND PLANE OF FIGURE 4, - SCHEMATIC OF A THIN WALLED COMPOSITE 

ISOTROPY. TUBE LONGITUDINALLY <0 « 0°) AND HELICALLY 

<E.G, 0 « 45°) REINFORCED, 
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figure 9 , - rorrm <e 1v c 22 ) m s hear <c 12 > strain VARI- 
ATION VERSUS FIBER ORIENTATION 0 ASSURING SAWLE 2 
RATER I Al PAAARETERS. 



60 80 100 


l 3 1 AND SHEAR it v ) STRAIN 
rATION 8 FOR SARPLE 2 
IE SHEAR STRESS STATE. 


ExIO' 



FIGURE 11, - NORMAL C£ n ,£ 22 ,E 33 > AND SHEAR <£ 12 > STRAIN 
VARIATION VERSUS FIBER ORIENTATION 0 FOR SAMPLE 2 
MATERIAL PARAMETERS UNDER OPPOSITE BIAXIAL STRESS STATE. 



FIGURE 12, - NORMAL AND RHew <£«> STRAIN 

VARIATION VERSUS FIBER ORIENTATION 0 FOR SAMPLE 2 
MATERIAL PARAMETERS UNDER TWO TO ONE BIAXIAL STRESS 
STATE, 



FIGURE 13, NORMAL (eL.E^.E^) AND SHEAR (e| 2 ) STRAIN 
VARIATION VERSUS FIBER ORIENTATION 0 FOR SAMPLE 2 
MATERIAL PARAMETERS FOR A CHANGE IN TEMPERATURE 
AT « 100 °F, 
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